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We investigate the spin structure of many-fermion systems 
with a spin-conserving two-body random interaction. We find 
a strong dominance of spin-0 ground states and considerable 
correlations between energies and wave functions of low-lying 
states with different spin, but no indication of pairing. The 
spectral densities exhibit spin-dependent shapes and widths, 
and depend on the relative strengths of the spin-0 and spin-1 
couplings in the two-body random matrix. The spin structure 
of low-lying states can largely be explained analytically. 



I. INTRODUCTION 

Random matrix theory has been widely used to model 
the behavior of complex quantum systems |l],|^] . Of par- 
ticular importance are two-body random matrix ensem- 
bles (TBRE), which were introduced into nuclear physics 
three decades ago and recently have received a con- 
siderable amount of renewed interest in the study of 
mesoscopic systems [PHlOj and nuclei juj-jlfl]. This is 
partly due to the observation of regular structures in the 
low-lying states of various nuclear models with random 
two-body interactions. As examples we mention the ob- 
served ground state dominance of J p = + states, the 
occurrence of pairing (energy) gaps in random shell mod- 
els fll2|| , and the existence of vibrational and rotational 
bands in random interacting boson models fll4]j . While 
these observations concern the low-lying levels of many- 
body systems, similar regularities also have been reported 
in the wave function structure. These include signatures 
of phonon collectivity relating states of different angular 
momentum 12 , strong pair-transfer amplitudes between 
ground states of nuclei with different mass numbers [fl7| , 
strong electromagnetic transitions in bands of random in- 
teracting boson models plM, deviations from random 
matrix theory in transition strength distributions for em- 
bedded Gaussian ensembles | jl3| , and regularities in wave 
functions of spin systems with random interactions [fl8|| . 
These observations indicate that some typical properties 
of interacting many-body systems are determined sim- 
ply by the presence of a rotationally invariant two-body 
force that allows for transitions between single-particle 
Fock space states. 

In this article we aim at an understanding of some of 
the reported regularities. We are particularly interested 
in the quantum numbers of the ground state and in corre- 
lations between energies and wave functions of low-lying 



states. For this purpose we consider a system of N spin-i 
fermions on M orbitals, subject to a spin-conserving ran- 
dom two-body interaction. This system is simpler than 
the nuclear models studied so far, and allows us to de- 
rive analytical results. It may be seen as a model for 
mesoscopic systems like quantum dots. 

This article is divided as follows. In the following Sec- 
tion we focus on the width and shape of spectral den- 
sities in sectors of definite spin for two models of spin 
systems with random interactions, and derive analytical 
results. Special emphasis is placed on understanding the 
spin structure of the spectrum close to the ground state. 
Our analysis extends and complements previous work in 
nuclei JlJ] and mesoscopic systems JhJ. In Section [II 
we investigate correlations between low-lying states, and 
consider pairing properties. Finally we give a summary. 



II. SPIN STRUCTURE 

A. The model 

We begin with a very simple model of random spin- 
conserving two-body interactions in a many-body sys- 
tem [HH}7): 

H = C J2 H aa ,AlA a ,+d HppA^Ap, 

a, a' spinO (3.(3' spinl 

(i) 

where a and a' each denotes a pair of two spin— i 
fermions in a J = state, while /3 and (3' each de- 
notes a pair of fermions in a J = 1 state. The 
two-body states can of course be labeled by the indi- 
vidual orbitals occupied by the two particles. Thus, 
for a system with M single-particle orbitals, we have 
M(M — l)/2 pairs with total spin J = 1, enumerated 
as Ap ( . = (aiia,j-\ + a^aj\) j where aj annihilates a 
single particle state, and the single-particle indices are 
1 < i < j < M. [Here we are enumerating two-body 
J = 1 states with J z = only; obviously there are equal 
numbers of J z — ±1 two-body states which by rotational 
symmetry must couple in precisely the same way.] Simi- 
larly, there are M(M — l)/2 pairs with total spin J = 0, 
having the form A aij — (a^a^f — a^aji)/\^2, and an ad- 
ditional M pairs where a single site is doubly occupied: 
A aii = aiia^. We will assume throughout that the two 
kinds of J = states interact equally strongly among 
themselves and with each other; this assumption does 
not substantially affect any of the results to be presented 
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in this paper. The effect of diagonal spin-spin interac- 
tions in a spin system with random wave functions has 
been investigated in ref. [ ^0| . 

Within each two-body spin sector, then, the couplings 
H aa i and Hpp* are assumed to be governed by a GOE 
random-matrix ensemble, 
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where (. . .) indicates an ensemble average. We note that 
the two matrices H aa ' and Hppi are each taken to be real 
and symmetric, consistent with time reversal invariance; 
the model can easily be generalized to the case where an 
external magnetic field is present, so that the matrices 
must be taken from a GUE rather than a GOE ensem- 
ble, or from an intermediate ensemble jl],^,|l4|]. Finally, 
we notice that in the model (Eq. [I]) , couplings in different 
two-body spin sectors are taken to be c ompletely inde- 
pendent of one another; in Section \1 J| we will discuss 
an alternative random-coupling model where this is not 
the case. In Section [II we will see what effects on spin 
spectra may be obtained by adding a one-body term to 
the two-body random Hamiltonian. Note that the pres- 
ence of a strong one-body term enhances the likelyhood 
of finding a spin-0 ground state, due to spin-degeneracy 
of the one-body spectrum. 



B. Preliminary results 

We may now consider numerically the spectra of 
N— fermion systems interacting via the random two-body 
interaction of the form (Eq. |l|) . Due to an approximate 
particle-hole symmetry of this system, maximum den- 
sity is obtained near half-filling, N — M , and we may 
restrict ourselves to the parameter range 2 < N < M. 
The preponderance of 5 = many-body ground states, 
which has been observed previously in shell model calcu- 
lations, is seen to be an extremely strong effect in this 
simple random interaction model. For N = 4 fermions 
on M = 4 sites, with an exclusively spin-0 interaction 
(Co = 1, C\ = Q), we find for an ensemble of 2000 systems 
that the ground state has total spin 5 = in 99.8% of 
the cases. For an exclusively spin-1 interaction (C\ = 1, 
Co = 0), the ground state is observed to have total spin 
5 = for 99.75% of the systems in the ensemble (a result 
which is statistically indistinguishable from the previous 
case). This is despite the fact that in this system 5 = 
many-body states comprise only about 56% of the total 
number of states. For N = 4 particles on M = 5 sites, 
all 2000 computed ground states were observed to have 
total spin 5 = 0, independent of whether a pure J = or 
J = 1 two-body coupling was used. Very similar results 
are obtained at least until a system size of M = 12, be- 
yond which point the size of the many-body basis makes 
the diagonalization of a large ensemble of Hamiltonians 



prohibitively expensive. For N = 6 particles on M = 6, 
7, and 8 sites, we once again observe an almost 100% 
prevalence of 5 = ground states. 

How may we understand these initially surprising find- 
ings within the context of a completely random interac- 
tion model? We begin by considering the shape of the 
many-body density of states ps{E) in sectors of total 
many-body spin 5. It is known that for N 3> 2 par- 
ticles with a two-body interaction, the moments of the 
spectral density move away from those associated with 
a semicircle law, and approach those corresponding to 
a Gaussian shape 20 21]. Then, the ensemble-averaged 
spectra (ps{E)}, being centered at E = because of the 
symmetry of the Hamiltonian, must be described sim- 
ply by their squared widths, Tt$H 2 (where Trg denotes 
throughout the sum over all basis states in the sector of 
total spin 5, divided by the number of basis states in 
that sector). Our naive expectation (assuming that the 
spectrum can be described by a Gaussian shape well into 
its tail) is that ground state behavior should be associ- 
ated with the relative widths of the different spin spectra, 
a preponderance of 5 = ground states arising from a 
large value of Tro-ff 2 as compared to widths in the other 
sectors [jlOlfll). Indeed, fixing the number of particles at 
N = 4, our numerical results show that for 4 < M < 10 
orbitals the relation Tr H 2 > TriiJ 2 > Tr 2 ff 2 is almost 
always observed (independent of the relative values of Co 
and Ci), consistent with 5 = ground state dominance. 
For M = 11 and 12, however, when taking an exclusively 
J = 1 coupling, this pattern is reversed, with the widths 
of the higher-spin spectra beginning to exceed those of 
the 5 = spectrum. This shows that the width analysis 
alone is not sufficient to explain numerical results for the 
ground state behavior. However, these findings do sug- 
gest than an interesting transition may be taking place 
as we approach the dilute regime, which we shall now 
proceed to understand analytically. 



C. Width analysis 

The many-body basis states of the system may be clas- 
sified by the total spin 5 and the number D of orbitals 
that are doubly occupied, < D < N/2. We easily see 
that in the dilute limit (M — > oo for fixed N), D = 
states will dominate: a typical 5 = or 5 = 1 basis 
state in this limit is given for even N by 
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where the plus (minus) sign produces an 5 = 1 (5 = 0) 
N— body state, the single-particle orbitals i, j, k\, l\, 
. . . , fcjv/2-i, In/2-i are all taken to be distinct, and |0) 
denotes the vacuum. However, at finite densities we must 
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consider basis states in subsectors of different D for any 
given value of S, as these might in general have spectra 
of different shape and width. It is a straightforward, 
though somewhat tedious, counting problem to compute 
the ensemble-averaged squared energy width (^\H 2 \^), 
for a typical basis state Y$>) of given quantum numbers 
S and £>, when acted on by the Hamiltonian (Eq. [I]). 
For each two-particle pair a 1 or (3' of the appropriate 
spin that may be annihilated, one simply has to count 
contributions from all pairs a or (3 that may be created 
in its place, while keeping track of terms that must be 
added up coherently. So, for example, for a pure J = 
coupling (Co = 1, C\ =0), one obtains 



Tr i? = 



iV 2 (2A/-iV) 2 2M 2 N + MN 2 - N 
64 + 

6MN + N 2 + 2N 2 



1G 



(jV/2) (iV/2) 
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, , 2M 2 - 2MN + N 2 
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In Eq. [|, terms are grouped and ordered by decreasing 
power in the number of orbitals M assuming finite den- 
sity (TV ~ M). In particular, we notice that the leading 
0(M 4 ) term is invariant under particle- hole symmetry, 
N — ► 2M — N. This symmetry is broken at 0(1/ M) be- 
cause of new terms that arise in the Hamiltonian (Eq. [I]) 
when one commutes annihilation operators A and cre- 
ation operators corresponding to pairs having at least 
one orbital in common. We notice also that the effect 
of double degeneracy (D > 0) is initially to increase the 
width slightly, and then restore it to the original value 
when D reaches its maximal value of N/2. However, the 
effect of double degeneracy appears only at 0(1/M 2 ), 
which we will see below to be small compared with the 
level spacing near the ground state, in the many-body 
limit. Thus, even for a non-dilute system (N ~ M), we 
are justified in studying the spectra of D = states and 
taking their widths as proxies for the entire spectrum of 
that spin, even for the purpose of determining behavior 
near the edge of the spe ctrum (where the level spacings 
are largest). In Section HI A we will consider the pos- 
sible relevance of D > contributions to wave function 
structure near the ground state. We find there that ba- 
sis states with high (or low) double occupancy have al- 
most equal overlaps with eigenstates of all energies, and 
have no special tendency to contribute to the low-lying 
part of the spectrum. Thus it appears that, to leading 
order, double occupancy effects are not important to un- 
derstanding the spin structure of the spectrum for a pure 
two-body interaction. 



D. The four-body system 

The width expression of Eq. ^ is rather unwieldy, as it 
consists of 14 monomial terms, and so we begin our de- 
tailed analysis with the special case of N = 4 particles; as 
we shall see below, all the interesting qualitative behavior 
of the many-body system is already present in this sim- 
ple case. In the S — sector we have 2( i f) independent 
basis states of type D = 0, having the form 
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(where i, j, k, and I all represent different orbitals), 



M( M 2 x ) basis states of type D 
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2. Similarly, in the S = 1 sector we have 
l j independent basis states of type D = 0, of the form 



- (at i0 } T + ot T o}J (a^ - 4 T 4j |0) , 



(0) 



and M( M 2 : ) basis states of type D = 1. Finally, in the 
5 = 2 sector we only have the (^) basis states of type 
D = 0. As the system becomes dilute (M — > 00), type 
D = states always dominate, and the total numbers of 
states in the spin 0, 1, and 2 sectors approach the ratio 
2:3:1. [As we mentioned above, even for finite M, the 
inclusion of D > states does not significantly affect the 
behavior of the spin spectra, leading to corrections that 
are small compared with the level spacing.] 

What then causes the preponderance of S = ground 
states in this system? Using D = basis states, one 
obtains the following expressions for the widths (of course 
the Cq part of these results is merely a special case of 
Eq. |): 



TrnH 2 = C 2 
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Trii? 2 = Co (M 2 - M) + C\ (2M 2 - AM + 4) 



(T) 



The leading 0(M 2 ) terms in the above expressions are 
easiest to understand. We simply need to take the num- 
ber of spin- J pairs in a four-body state of total spin S 
(Eq. |^ or Eq. ^|), and multiply by the number of orbital 
pairs to which the two particles may jump (w M 2 /2). 
So, for example, in a state of total spin S — (Eq. ||), 
the ij and kl pairs are in pure J = combinations, while 
the remaining 4 pairs are uncorrelated, i.e. 1/4 of them 
have J = and 3/4 have J = 1. This leads to a total 
of 3 pairs each of J = and J = L explaining the two 
|M 2 terms in the first line of Eq. |^. In a state of total 
spin S = 1 (Eq. |J) , the only difference is that the ij pair 
is now in a J = 1 combination, so there are only 2 pairs 
with J = and 4 pairs with J — I, leading to the CqM 2 
and 2C 2 M 2 terms in the second line of Eq. [f| We see 
already from this very simple argument that, due to a rel- 
atively larger number of J = pairs in a spin-0 state and 
a relatively larger number of J — 1 pairs in a spin-1 state, 
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in the M — » oo limit we always have TtqH 2 > TiiH 2 if 
C > Ci and Tr ff 2 < TriiJ 2 if C < C x . The ratio of 
the widths in the dilute limit approaches a constant, 



Tr iJ 2 _ 3C 2 + 3C 2 
Tni? 2 ~ 2C 2 + 4C 2 ' 



(8) 



while the ratio of the many-body level spacing near the 
ground state, p^^Egs), to the width and thus to the 
width difference is becoming small. Because the spectral 
shapes po,i(E) approach a Gaussian form in the M — *■ oo 
limit, we expect this width difference to dominate the 
ground-state behavior, and to see 100% 5 = ground 
states for a strong Co coupling, and 0% for a stronger C\ 
coupling. 



or S = 1 ground states (i.e. the ratio of S = to S = 1 
ground states should approach the 3 : 2 ratio valid for 
the total density of states). We note, however, that in 
our numerical simulations, the ground state always comes 
out to have spin 5 = 0, even for the highest values of M 
(M = 12) that we can achieve, where the width diffe rence 
is very small. Indeed, as we saw above in Section [ID, 
5 = ground states can dominate even when the width 
difference changes sign, as happens for a pure spin-1 two- 
body coupling (Co = 0). A major reason for this we 
believe to be strong deviations from the Gaussian shape 
of a four-particle spectrum (though correlations between 
Po(E) and pi(E) may also be relevant; w e will discuss 
these correlations in more detail in Section IIIB). 



E. The Co = Ci special case 

The M — > oo behavior for Co = C% is more subtle 
and requires an analysis of the subleading 0(M) terms in 
Eq. ^[ To compute these we need to count more carefully 
the number of places a or (3 to which the two interacting 
particles may jump, given that two of the orbitals remain 
occupied by the two 'bystanders' during the interaction. 
In the case of a J = 1 pair, we also need to treat cor- 
rectly the interference between J z — 0, ±1 terms. The 
end result is that overall, after adding together J = 
and J = 1 couplings, the 5=1 width is suppressed at 
O(M) compared with the 5 = width, as a simple phys- 
ical argument shows. In Eq. H, the jk pair (or any one 
of the other 3 analogous pairs) may interact either in a 
J = or J = 1 combination, with the final state leaving 
j doubled up with i and with k ending up on any of the 
remaining M — 1 sites. For an initial state of the form 
Eq. |(| this process is not allowed (since the ij pair start 
out in a spin-symmetric combination and cannot end up 
occupying the same orbital), so TriiJ 2 is reduced at the 
O(M) level. Thus, in the large-M limit for Co = C%, we 
obtain 



Trpff 2 
Trii? 2 



1 + 0(l/M) . 



(9) 



By itself, however, this width difference should not 
be sufficient to determine the ground state behavior. 
That is because for a Gaussian spectral density with d 
states the location of the ground state scales as E gs ~ 
VTri? 2 VTogd, while the level spacing near the ground 
state goes as p~ 1 (E gs ) ~ VTrH 2 /y/logd. So the ratio of 
the expected difference (gap) between the lowest 5 = 
and 5=1 states to the level spacing scales as 



log d 4 log M 
~M M 



0. 



(10) 



and thus in a model of random uncorrelated fluctuations 
under the two Gaussian spectral envelopes po(E) and 
Pi{E), we would expect no preference for either 5 = 



F. Spectral shape effects 

We recall that the approach of an N— body spectrum 
to a Gaussian shape given only two-body interactions de- 
pends crucially on the assumption that repeated applica- 
tions of the Hamiltonian act on different pairs of particles 
and therefore may be taken to commute with each other. 
So, for example, 
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where a, a', a", and a"' label two-body states (either 
J = or J = 1), as in Eq. ^, and are summed over. 
Notice that the third term contributes in general only 
if the two transitions involve distinct pairs of interact- 
ing particles. Otherwise (e.g., for a generic sparse matrix 
with no two-body interaction structure) , we would obtain 
(TiH 4 ) — 2(Tr_ff 2 ) 2 , the relation appropriate to a semi- 
circular spectral shape rather than to a Gaussian [^o] ^2| . 
So the extent to which the spectrum deviates from a 
Gaussian shape towards a semicircle may be estimated 
by considering the probability that successive two-body 
interactions do not involve distinct pairs of particles. 

Consider for example a pure J = 1 two-body coupling. 
In an 5 = four-body state (Eq. ||), there are a total 
of three J = 1 pairs, as we have already seen above. 
After two of the particles interact in a two-body J = 1 
combination, the remaining two particles are guaranteed 
to be also in a J = 1 state. So there is a 1/3 probability 
that the next interaction will involve the two particles 
that have not yet interacted, and thus (Tro-ff 4 ) w (2 + 
■|)(Tro-ff 2 ) 2 . On the other hand, if the total many-body 
state is either 5=1 (Eq. ||) or 5 = 2, similar counting 
arguments lead to the conclusion that there is only a 
1/6 probability that the second interacting J = 1 pair 
will be disjoint from the first; thus (Tr^ai? 4 ) ~ (2 + 
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f )(Tri ) 2-ff 2 ) 2 . This behavior for the spectral shapes has 
been numerically confirmed (e.g., for M = 10 orbitals, 

/^gsp = 2.42, 2.18, and 2.15, for total many-body spin 
5 = 0,1, and 2, respectively), and is similar for the J = 
two-body interaction. So in particular for the Cq = C\ 
case of equal couplings, although the widths of the 5 = 
and 5=1 spectra approach each other in the M — > oo 
limit (Eq. |), the fourth moments remain in a finite ratio 
in this limit and lead to po(E) having a substantially 
longer tail than pi{E). Therefore, despite the almost 
equal widths, ground states in a Co = C\ system will be 
predominantly 5 = 0, even in the dilute M — + oo limit. 



G. The general many-body problem 

The above analysis of the N = 4 particle system gen- 
eralizes easily to the real iV-body problem that we are 
interested in. We begin, as before, by addressing the be- 
havior of the spectral widths in the 5 = and 5=1 spin 
sectors. Typical basis states in these two sectors were al- 
ready given above in Eq. [| Now consider a pure J = 
coupling. From Eq. || we get 



Trpff 2 
Triff 2 



l + O 
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iV 2 



(12) 



in the dilute M — > oo limit. This can be understood as 
follows. The only difference scaling as M 2 in the spectral 
width is that the ij pair in Eq. |] is in a J — combination 
for a many body-state of total spin 5 = 0, but not if the 
many-body state has total spin 5=1. Thus, the 5 = 
state has one more J — pair that can be acted on by the 
Hamiltonian, out of a total of w N 2 /8 pairs with J = 0. 
Similarly, for a pure J — 1 coupling we have one extra 
pair that can be acted on in an 5 = 1 many-body state 
(out of a total of « 3iV 2 /8 pairs with J = 1), so finally 
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If we keep the number of particles N fixed and take the 
system size M — ► oo, then the many-body level spac- 
ing of course goes to zero, while the width ratio remains 
finite for Co ^ C\. Once again we see, this time for gen- 
eral values of TV, that assuming the spectral structure is 
dominated by the spectral width, we expect 100% (0%) 
5 = ground states if Co > Ci (Ci > Co), in the dilute 
M — » oo limit. 

However, we saw in our discussion above of the N = 
4 special case that spectral cumulants beyond TtsH 2 
need also to be included in a full analysis. The ratio 
TrsH 4 /(TrsH 2 ) 2 always favors a longer tail for 5 = 0, 
so we only need to consider the case of a J = 1 coupling, 
where the width effect and the spectral shape effect are 
competing with one another. Referring back to Eq. 
we note that if particles on orbitals i and k z interact 



via a two-body J = 1 interaction, then the two parti- 
cles which remain on sites j and l z are also in a J = 1 
state and can subsequently interact if and only if the to- 
tal system is in an 5 = state (so that the four particles 
i, j, k Zl l z are in an 5 = state). Because the j, l z 
pair is distinct from the i, k z pair by construction, this 
causes it to be more likely that the next J = 1 pair of 
interacting particles will be distinct from i, k z if S = 0, 
and contributes towards an enhancement of the shape ra- 
tio Tr 0j ff 4 /(Tr # 2 ) 2 as compared with Tri J ff 4 /(Tri J ff 2 ) 2 . 
However, in a many-body system this effect is suppressed 
by 0(1/7V 3 ) (because there are 0(N 4 ) ways to choose 
two pairs, and only O(N) of them are involved in this 
effect). So for Ci > C the 0(\/N 2 ) difference in widths 
between p\ (E) and po (E) (Eq. [j^) is expected to exceed 
the difference in spectral shapes, leading to a dominance 
of 5 > states near the edge of the spectrum, in the di- 
lute many-body limit (which we define as M — > oo first, 
followed by N — > oo). 

For Co > Cl, on the other hand, both the difference in 
widths and difference in spectral shapes favors an 5 = 
ground state, and this behavior is expected to persist into 
the dilute many-body limit. For Co = Ci, the relative 
width difference vanishes (Eq. [l3]) in this limit, but the 
difference in spectral shapes now takes over, again leading 
to 5 = dominance near the edge of the spectrum. 



H. Transition at finite densities 

We have seen that for finite (large enough) number of 
particles N, 5 = ground state dominance is expected 
in the dilute M — > oo limit, as long as the spin-zero 
two-body coupling is at least as large as the spin-one 
coupling, Co > Ci. In the contrary case, the situation 
is reversed, with higher-spin states dominating behavior 
near the edge of the spectrum. What happens, though, 
if the system is not very dilute? For a finite number of 
single-particle orbitals M, we must take into account the 
subleading behavior in M of the spectral widths. An 
argument very similar to the one leading to Eq. but 
generalized to large N and to an arbitrary ratio of cou- 
plings Ci / Co shows that the effect discussed there always 
favors the 5 = width, at 0(1/NM), independent of 
which coupling is stronger. So for Co > Ci nothing par- 
ticularly interesting happens at finite densities; the pref- 
erence for 5 = states near the edge is only enhanced by 
finite density effects. This preference, however, becomes 
small (compared to the level spacing) in the many-body 
limit. 

On the other hand, for a stronger spin-1 coupling 
(Ci > C ), we expect the 0(l/iV 2 ) effect favoring 5 = 
ground states and the 0(1/NM) effect favoring 5 > 
ground states to balance each other, leading to a transi- 
tion in the width ratio at finite densities. In fact a closer 
analysis of the 4— body system (Eq. 0), setting equal the 
terms proportional to C\ in the two expressions, shows 
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that for a pure J = 1 coupling we have a transition just 
below M = 11, with a larger 5 = width for M < 11 and 
larger 5 > widths for M > 11. This is consistent with 
our numerical observations above. For comparable cou- 
pling strengths, the transition of course will take place 
at lower densities; so, for example for C\ — Cq(1 + e), 
the transition in the width ratio occurs near M = 10/e. 
Although we are unable to observe this transition nu- 
merically for larger numbers of particles, the analytical 
arguments suggest that one should take place in general 
at densities N/M ~ e for a small coupling difference e and 
at an O(l) finite density for a pure J = 1 interaction. 



I. Numerical results 

In Table [| we show some results of a numerical calcu- 
lation for N = 4 particles distributed over 5 < M < 12 
orbitals, with a pure J = 1 two-body coupling (C\ = 1, 
Co = 0). The first row of the Table shows the ratio of 
widths \/Tr sH 2 in sectors of total many-body spin 5 = 
and 5 = 1. We see that Tt H 2 > Ty x H 2 holds at low M 
(high density), while in the dilute regime (M > 12) the 
relation is reversed. The crossover occurs near M = 11 
for N = 4, in good agree ment with our analytical expec- 
tations (see Section [IH). The same qualitative behavior 
is predicted and observed for the higher spin sectors; we 
see from the second row of the Table that the crossover 
between 5=1 and 5 = 2 dominance occurs at M = 9. 

We proceed to compare the energies of the lowest-lying 
states in the different spin sectors, and find a trend to- 
wards lower-lying states of higher spin as M increases, en- 
tirely consistent with the results for the spectral widths. 
However, by comparing rows 1 and 3 (or 2 and 4) of the 
Table, we see that for any given system size M, low-lying 
states with small total spin are preferred more than would 
be expected based simply on the width ratio. Thus, for 
M = 11, where the spin-0 and spin-1 widths are nearly 
equal, the ground state of the system is still nearly al- 
ways 5 = 0, and on average the ground state is observed 
to have energy 12% lower than that of the lowest 5=1 
state. This is due largely to th e spectral shape effects dis- 
cussed above in Section II F , which favor ground states 
of lower total spin given equal spectral widths in the dif- 
ferent spin sectors. 

When comparing 5 = and 5 = 1 sectors in 
particular, 



we 



recall from Section II F that 
TrgiJ 4 / (Tr sH 2 ) 2 is significanlty larger for 5 



the ratio 
= 0, for 

any value of M. For N > 4 (see Section II G), spectral 
shape effects are reduced, but unfortunately we are un- 
able to probe the dilute M ^> N regime numerically for 
N > 6. Instead we stay at N = 4 particles and compare 
the 5=1 and 5 = 2 sectors, which as we saw previously 
have nearly equal values of the ratio Trgff 4 / (TisH 2 ) . 
Indeed, by M = 12 the mean values of the lowest-lying 
5 = 1 and 5 = 2 states become almost equal, and al- 
ready by M = 9 the lowest-lying 5 = 2 state appears be- 



fore the lowest-lying 5=1 state in a significant fraction 
of systems in the ensemble. The last row of the Table 
supports the analytical arguments presented above for 
a transition between lower-spin and higher-spin ground 
states as a function of the particle density. 



J. An alternative model 

As we have mentioned previously, the model (Eq. |l|) 
assumes that there is no correlation between two-body 
interactions in J = and J = 1 two-body states. We 
may instead consider an alternative model, motivated by 
Ref. pcfl and much subsequent work, in which the inter- 
action Hamiltonian is separated into a spin-independent 
and a spin-dependent part: 



H = J2 J2 [CvV ljkl + C w W ijkl J k ■ J, 

ijkl JiJjJ k Ji 



X a iJi a jJ. a kJ h a UiS(Ji + J J ~ Jk - Jl) 



(14) 



where the indices i, j, k, and I label individual orbitals, 
Ji, Jj, J k , and Ji = ±s are the z— components of spin 
for each of the four fermions, and total spin conservation 
has been enforced. The matrix elements Vijki and Wij k i 
are taken to be real Gaussian random and independent 
variables, with the constraints 



Vuik — Vikjl — Vwd — Vi 



Ijki 



•'ijkl 



(15) 



(and similarly for the W%jki) imposed by the reality of 
wave functions and the symmetry of the two-body inter- 
action. The constants Cy and Cw determine, of course, 
the relative strengths of the spin-independent and spin- 
coupled interactions. 

The behavior near the ground state for the model of 
Eq. ^ can be understood in terms very similar to those 
we used for analyzing the original model of Eq. [I]. We 
assume initially that spin-spin interactions cannot be ne- 
glected, i.e. Cw 7^ 0. We first note that the spin- 
spin interaction J k ■ J t = (l/2)((J fe + J/) 2 — J 2 — J 2 ) — 
(l/2)((Jfc + Ji) 2 — 3/2) is three times stronger when par- 
ticles k and I are in a J = two-body state as com- 
pared with a J = 1 two-body state. Now we have seen 
above that an 5 = TV— particle many-body state has 
0(1/N 2 ) more J = pairs than an 5 = 1 many body 
state: this leads to a wider po(E) envelope and 5 = 
ground state dominance, just as for the Co > C\ situa- 
tion considered in the previous model. Thus, a spin-spin 
interaction (of either sign) is seen immediately to lead 
to 5 = ground state dominance (100% dominance in 
the M — > oo many-orbital limit). At finite densities, 
N/M ~ const, we must once again include corrections to 
the wid th d ifference, of order 1/NM, as in Eq. ^ and in 
Section [IH above. We find, however, that the subleading 



correction, associated with the impossibility of doubling 
up two particles in a J — 1 two-body state, always favors 



G 



the S = many-body state. This is a model-independent 
property resulting simply from the counting of two-body 
J = vs. J = 1 pairs in a many-body state. Therefore 
the 0{1/N 2 ) and 0(1/NM) effects contribute with the 
same sign, and even at finite density we expect an excess 
of S = ground states. Finally for sma ll systems, shape 
effects (discussed above in Section [IF) also contribute, 
again favoring the S = ground state. 

We may also consider the scenario of a spin-free inter- 
action (Cw = 0), which because of the equal strength 
with which it couples S — and 5 = 1 pairs may be 
compared with the Co = C\ special case discussed above 
for the model of Eq. |l]. As in that case, the width dif- 
ference is now expected to vanish in the M — > oo dilute 
limit; furthermore as before the width difference eventu- 
ally becomes small compared with the level spacing in 
the system. However, at moderate values of the parti- 
cle number N, spectral shape effects are expected to be 
very important, leading to S — ground state excess at 
arbitrarily large values of M. 



III. CORRELATIONS AND PAIRING 

In realistic many-body systems there are strong corre- 
lations between ground states of different quantum num- 
bers, and the pairing correlation is of particular impor- 
tance. This correlation has also been observed in various 
studies of nuclear models with random two-body inter- 
actions [ p^ , pT|Jr^| . In this Section we investigate corre- 
lations between energies and wave functions of low lying 
states with different quantum numbers. 



A. Pairing 

Pairing is a well known phenomenon in many-body 
fcrmion systems. Its spectral signature is an energy gap 
that separates the ground state from the excited states. 
In the presence of pairing, the ground state is a conden- 
sate of pairs of fermions whose quantum numbers are 
related to each other by time reversal symmetry. In the 
case of the spin systems considered in this work, pair- 
ing refers to a ground state of the form JJ . b^b'^ |0), 
where &t is a quasi-particle operator that is related to 
the single-particle operators a\ by a unitary transforma- 
tion, as bj = 2i=i Ujidi. 

Let us first consider pairing in the single-particle ba- 
sis generated by the a\. The results on wave function 
structure in random spin systems reported in Ref. JlSj ] 
show that states at high spectral densities have a num- 
ber of principal components that agrees with standard 
random matrix theory while low-lying states have a con- 
siderably smaller number of principal components (thus 
indicating a simpler structure and some degree of regular- 
ity I23J24]). However, the structure even of the low- lying 



states was found to be much more complex than a sim- 
ple Slater determinant. Based on that previous analysis, 
and on the analytical results obtained above in Eq. ^ 
showing very weak dependence of spectral width on dou- 
ble occupancy number, we do not expect to find pairing 
here in the single-particle basis. To compare these ar- 
guments with numerical data we measure the number of 
doubly occupied single-particle orbitals in a given eigen- 
state \E, S) with energy E and total spin S by comput- 
ing P E<S = (E, S\J2j a li a ]f a n a ji\ E i S) ( com pare the 
related quantity D for basis states in Section [I C ) . For an 
ensemble of 100 random Hamiltonians with N = 6 spin- 
| fermions on M = 6 sites, Co = 1, and C\ = 0, Pe,s 
shows fluctuations around {Pe,o} ~ 1.71, (Pe,i) ~ 1-43, 
{Pe,2) ~ 0.86 and does not exhibit any significant energy 
dependence. Thus, no pairing is present in the single- 
particle basis a\. The numerical results agree well with 
the analytical expressions Pe.o = 12/7, Pes = 270/189, 
and Pe,2 — 6/7. 

Pairing in a general quasi-particle basis can most con- 
veniently be investigated by examining its spectral sig- 
nature. The existence of a pairing gap is closely re- 
lated to the approximate degeneracy of the first excited 
spin-0 level E^~ ^ with the lowest-energy spin-one state 

(5—1) 

E\ . Within the two different models considered in 
this work, we do not observe any degeneracy between 
E^ -0 ^ and e[ S_1 \ and find the ratio of splitting to gap 
as 



Aee 







=1) 




min 


[e^\e[ s ^) 





2.8 



(16) 



for an ensemble of 100 Hamiltonians with Co = 1, Cl =0, 
and M = N = 6. Similar results hold for other param- 
eter values, and we never see A « 1, which would be 
indicative of a pairing structure. Thus, the dominance 
of spin-0 ground states cannot be caused by pairing cor- 
relations within our model. Note that this result differs 
from the observations made for the nuclear shell model 
with random interactions jl2|,[l7],[l4]] . As an illustration 
we show in Fig. I the spin dependent spectra for a sys- 
tem of six fermions on six single-particle orbitals. The 
random two-body interactions only couple pairs of spin 
J = 0, i.e. we set C\ = 0. The ground state and the 
low energy states have spin S = 0. The absence of low- 
lying spin-degenerate states clearly indicates the absence 
of pairing correlations. 

To study further the subject of pairing, we add a one- 
body part of the form 



M 



(17) 



to the Hamiltonian (Eq. [j]). Fig. ^ shows the pairing ob- 
servable A as a function of c = C/(NMGq) for C\ = 
and N = M = 6. For c> lwe find A <C 1, indicating 
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strong pairing. In this regime the two-body interaction is 
a small perturbation, and we obtain pairing in the basis 
of single particle states. Note that this behavior is not 
caused by the harmonic form of our one-body potential 
but rather is determined by the applicability of a single- 
particle picture. For decreasing values of c one leaves 
this perturbative regime, and the increase of A shows 
that pairing is suppressed. Note that A > 1 at suffi- 
ciently small values of c. This indicates that there are 
several spin-0 levels below the spin-1 ground state, as in 
Fig. ^] above. This finding is consistent with the results of 
Section |l| and hints at the dominance of spin-0 low-lying 
states. As a check we also replaced the spin-1 ground 
state energy by the highest spin-1 energy (which should 
not be at all correlated with the spin-0 ground state en- 
ergy), while correcting for a nonzero mean of the total 
spectral density, i.e. replacing £^ 5-1 ^ — > TrH — , 
and recomputed A. For the system with a pure two-body 
force we found no difference when compared to the pre- 
vious results, confirming that the value of A at c = is 
consistent with random spectral fluctuations. 



B. Level and energy correlations 

The absence of pairing in the presence of a strong two- 
body force is peculiar to our model and is not shared 
by nuclear shell models with random interactions. It is 
therefore interesting to look for different correlations be- 
tween low lying states of spin-0 and spin-1. Such correla- 
tions have been observed previously in Refs. |12 17|. We 



therefore look at the generalized collectivity defined as 



\(S = 1\B\S = 0)\ 2 
1 (S = Q\B*B\S = QY 



(18) 



where B = Y,i<j c ij B i] an d 

Bij = (a- ; £i] T + at T at i )(a i4 ,aj T - a lT a iX ) (19) 

is a spin-flip operator that transforms a pair of spin zero 
into a pair of spin one. Following Refs. jl2|,|l^| we choose 
the coefficients as 



Cij = (S=l\Bij\S = 0). 



(20) 



This is more convenient than searching for those coeffi- 
cients that maximize /. To be definite, we take the 
ground states of spin zero and spin one, i.e. \S = 0) = 
\e[ S=q) ) and \S = 1) = \E[ S=1) ), and compare with 
the corresponding value that is obtained by substitut- 
ing \S = 1) — > |.Emax )• The states corresponding to the 
highest and lowest energies in a sector of definite spin are 
assumed to be uncorrelated. The results are presented in 
Table |l| for N = 4 and different values of M, once again 
setting C\ = 0. Obviously, there is a strong correlation 
between the ground states of different spin. It is instruc- 
tive to let \S = 1) run through the n lowest lying spin- 



\S = 0) = \e[ S ~^) fixed. Figure^ shows that the corre- 
lation is strongest for the ground states and decreases for 
the higher lying states. This finding is consistent with the 
recent observation Jl8| that low-lying states are weakly 
localized in Fock space, while high-lying states are delo- 
calized and agree with random matrix predictions. 

It is also interesting to consider energy correlations be- 
tween ground states of spin zero and spin one. The ap- 
propriate correlator is 



G 



((^= o) ) k =i) )> 



(s4 



(S=0) 



5E, 



(S=l) 



1/2 



(21) 



where SE^ = E^ — (^Eq S ^ and (. . .) denotes an en- 
semble average. We find G = 0.93 for N = 4 particles on 
M = 9 sites (with C\ =0), indicating a strong correla- 
tion. The observable G, is however, somewhat mislead- 
ing, since a significant part of the fluctuations 5Eq from 
one realization of the ensemble to another arises simply 
from the fluctuation of the centroid and width of the over- 
all spectral envelope, TtH and TvH" . The centroid and 
width for a given system are almost identical in sectors 
of different total many-body spin. To understand fluctu- 
ations under the spectral envelopes, we should therefore 
remove this trivial effect by considering the correlation 
G of E { S=0) and E { S=1 \ where 



E 



E - TrH 



TrH 2 - (TtHY 



(22) 



1 states \S = 1) = \E{ 



(S=l)\ 



\E, 



(S=l)\ 



while keeping 



We then obtain G = 0.76 for the same set of parameters. 
The value of G increases as one goes to larger systems; 
e.g. for N — 4 particles, G increases monotonically from 
0.20 to 0.76 as M goes from 4 to 9. These observations 
are consistent with the existence of correlations between 
the wave functions of the corresponding states. 

Note that the observation of strong correlations be- 
tween low-lying states of spin-0 and spin-1 is in qualita- 
tive agreement with results in shell models and interact- 
ing boson models. However, in contrast with the findings 
in these nuclear models, we do not find any indications 
of a pairing correlation in the spin systems considered in 
this work. One can only speculate that this is due to 
the additional isospin degrees of freedom and the more 
sophisticated angular momentum coupling scheme used 
in the nuclear models. 



IV. SUMMARY 

We have investigated the spin structure of low lying 
states in fermion systems with random two-body inter- 
actions. Numerical calculations show that the ground 
state has spin zero with almost 100% probability. For 
a dilute system we derive analytical expressions for the 
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shapes and widths of the spectral densities at different 
spin. These explain our numerical results. For systems 
with a spin-1 coupling that is stronger than the spin-0 
two-body coupling we predict a transition from spin-0 
ground states to ground states of higher spin. While it is 
difficult to access this regime numerically, we do observe 
one of its precursors - namely a spin-2 ground state that 
is lower in energy than the spin-1 ground state. 

Furthermore, we have shown that there are strong cor- 
relations between wave functions and energies of low- 
lying states of different spin quantum number. The wave 
function correlations are manifest in expectation values 
of spin transition operators and indicate some degree of 
collective behavior. The existence of energy correlations 
confirm the picture about the dominance of spin-0 ground 
states. These correlations decrease as one moves away 
from the ground state region. In the absence of a one- 
body force we find no indication for any pairing corre- 
lation between low lying states, and our model differs 
herein from the nuclear models with random interactions. 
Pairing correlations are only observed after a sufficiently 
strong one-body potential is added to the random Hamil- 
tonian. 
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M 


5 


6 


7 


8 


9 


10 


11 


12 


a (6'=0) /(T (5=l) 


1.20 


1.13 


1.09 


1.06 


1.03 


1.01 


1.00 


0.99 


a (6'=l) /(T (5=2) 


1.21 


1.12 


1.07 


1.03 


0.99 


0.97 


0.95 


0.94 


£f= UJ /E<f =1> 


1.40 


1.27 


1.21 


1.17 


1.14 


1.11 


1.12 


1.09 


E V*=V/ E V>=*) 


1.89 


1.40 


1.21 


1.14 


1.07 


1.03 


1.02 


1.01 




0% 


0% 


0% 


0% 


0% 


0% 


0% 


0% 


E (S=2) < E (S=1) 


0% 


0% 


0% 


< 1% 


10% 


20% 


20% 


70% 



TABLE I. Relative spectral widths cr (s) = VTrsif 2 and 
energies of the lowest-lying states E{ in sectors of total 
many-body spin S — 0, 1, and 2, for N = 4 particles on 
M sites, with a pure J = 1 two-body coupling. The last 
row shows the fraction of systems in the ensemble in which 
the lowest-lying S = 2 state is lower in energy than the low- 
est-lying S — 1 state. 



M 


4 


5 


6 


7 


8 


9 


f 

,/random 


0.45 
0.11 


0.37 
0.05 


0.33 
0.03 


0.29 
0.02 


0.24 
0.01 


0.26 
0.01 



TABLE II. Wave function correlation / between ground 
states of total spin zero and spin one in a system of N — 4 
spin-i fermions on M orbitals. /random measures the corre- 
lation between the (presumably uncorrelated) spin-0 ground 
state and the highest state of spin one (see text for details). 
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FIG. 1. Spectra for a system of six fermions on six orbitals 
with S z — as a function of total spin S. 
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FIG. 3. Correlation / (see text) between the spin ground 
state and the n th state of spin 1. 
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FIG. 2. Pairing in a many-body system with a random 
two-body interaction and a harmonic one-body potential. A 
measures the ratio of the splitting to the gap, while C/Co 
measures the ratio of the one-body and the two-body cou- 
plings (see text for details). 
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